Draft version January 14, 2013 

Preprint typeset using IAT^X style emulatcapj v. 04/03/99 



THE FORMATION OF THE ECCENTRIC-ORBIT MILLISECOND PULSAR J1903- 
THE ORIGIN OF SINGLE MILLISECOND PULSARS 



-0327 AND 



O 

(N 



43 

Q-i! 
6 



> 
in 

m 

cn 
o 



S. PORTEGIES ZWART 1 , E.P.J. VAN DEN HEUVEL 2 , J. VAN LEEUWEN 3 , AND G. NELEMANS 4 

1 Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The 

Netherlands 

2 Astronomical Institute 'Anton Pannekoek', Science Park 904 1098 XH Amsterdam, 

The Netherlands 

3 Stichting ASTRON, PO Box 2, 7990 AA Dwingeloo, The Netherlands 
4 Dept. of Astrophysics, Radboud University Nijmegen, Heyendaalseweg 135, NL-6525 

AJ Nijmegen, the Netherlands 

Draft version January 14, 2013 
ABSTRACT 

The millisecond pulsar J1903+0327 is accompanied by an ordinary G-dwarf star in an unusually wide 
(P or b — 95.2 days) and eccentric (e ~ 0.44) orbit. The standard model for producing MSPs fails to 
explain the orbital characteristics of this extraordinary binary, and alternative binary models are unable 
to explain the observables. We present a triple-star model for producing MSPs in relatively wide eccentric 
binaries with a normal (main-sequence) stellar companion. We start from a stable triple system consisting 
of a Low-Mass X-ray Binary (LMXB) with an orbital period of at least 1 day, accompanied by a G-dwarf 
in a wide and possibly eccentric orbit. Variations in the initial conditions naturally provide a satisfactory 
explanation for the unexplained triple component in the eclipsing soft X-ray transient 4U 2129+47 or 
the cataclysmic variable EC 19314-5915. The best explanation for J1903 + 0327 however, results from 
the expansion of the orbit of the LMXB, driven by the mass transfer from the evolving donor star to 
its neutron star companion, which causes the triple eventually to becomes dynamically unstable. Using 
numerical computations we show that, depending on the precise system configuration at the moment 
the triple becomes dynamically unstable, the ejection of each of the three components is possible. If the 
donor star of the LMXB is ejected, a system resembling J1903 + 0327 will result. If the neutron star 
is ejected, a single MSP results. This model therefore also provides a straightforward mechanism for 
forming single MSP in the Galactic disk. We conclude that the Galaxy contains some 30-300 binaries 
with characteristics similar to J1903 + 0327 and about an order of magnitude fewer single millisecond 
pulsars produced with the proposed triple scenario. 



1. INTRODUCTION 

The classic channel for the formation of a millisecond 
pulsar (MSP) requires a close binary with an extreme mass 
ratio ( £ 2/10). This binary survives a common-envelope 
evolution and the subsequent supernova explosion of the 
primary star. The neutron star, formed in the supernova, 
can subsequently be spun up to a millisecond pulsar (MSP) 
in a phase of mass transfer from the ^ 2 M Q Roche- 
lobe filling companion star, which in the process is slowly 
stripped from its envelope. During this phase the binary is 
visible as a low-mass x-ray binary (LMXB), eventually re- 
sulting in a MSP that is accompanied by a low-mass white 
dwarf in a relatively wide, almost c ircula r orbit (e £ 10 -3 ) 
( Bha ttacharva fc van den Heuvell . [l991). In the Galaxy 
50 such systems are known, while about 15 MSP s in the 
Galaxy have no companion at all (|Lorime n. l2008h . 

The formation of the recently observed binary millisec- 
ond pulsar J1903 + 0327 cannot be reconciled with the 
above scenario. It's characteristics are too different: the 
companion star is a G-dwarf instead of a white dwarf, 
and the orbit is hig hly eccentric, e ~ . 44 ins tead of the 
expected < 10~ 3 (jCordes fc Chernofil 1 19981 ). In addi- 
tion, the average mass of the companion of known Galac- 
tic MSPs with pulse period < 10ms is 0.22 ± 0.17M©, 
whereas for J1903 + 0327 the companion mass is 1.03 M© 
iFreire et al] (|2010l ); and binary MSPs with a companion 



mass > 0.6 Mq tend to have a long pulse period (P) = 
62±76ms and short orbital periods (P rb = 5.2±4.8days), 
whereas J1903 + 0327 has an extremely short pulse pe- 
riod of 2.15 ms and an extraordinary long orbital period 
of 95 days. These discrepancies with respect to the ex- 
pected outcome of the standard scenario for producing a 
millisecond pulsar in a binary requires an exotic solution. 

We propose that J1903 + 0327 was born as a rather 
ordinary triple star of which the inner binary is the pro- 
genitor of a LMXB, and with an outer (tertiary) star that 
initially is less massive than the secondary so that the 
inner secondary evolved first. After a common-envelope 
phase and a supernova explosion, mass transfer in the in- 
ner binary leads to expansion of its orbit. Depending on 
the orbit of the outer star a dynamical instability ensues 
in which one of the three stars is ejected. Such an evo- 
lution results in a MSP binary with an outer companion 
in either a wide orbit (if the instability is avoided), or a 
MSP binary if the outer companion is ejected, or a normal 
binary plus a single MSP if the neutron star is ejected, or 
if the inner secondary is ejected, a MSP with a low-mass 
companion in an eccentric orbit, like J1903 + 0327. 

The scenario sounds exotic, and it is, but in §[3] and 
§0 we estimate the available parameter space and cal- 
culate that this model leads to a n acceptable birthr ate. 
(We notice that in a recent paper IFreire ei~al\ (|201dD in- 
dependently also suggest this model one of possible triple 
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star models for J1903 + 0327.) We will now dwell on the 
details of our scenario in §[5] and quantify the results by 
performing simulations of triple star systems in §|4] We 
discuss the shortcomings of earlier proposed scenarios for 
the formation of J1903 + 0327 in §[B] and summarize our 
conclusions in §[7] 

2. THE TRIPLE SCENARIO FOR FORMING J1903 + 0327 

We propose that the binary J1903 + 0327 was born as 
a rather ordinary triple star. We envision a configura- 
tion at birth where a relatively massive 9-12 Mq primary 
star was orbited by a secondary star of ~ 0.8 to 2.0 M Q 
in a relatively close (~ 200 R©) orbit, and a tertiary star 
which is less massive than the secondary and has a rather 
wide ^ 560 R© orbit around the inner binary. Adopting a 
tertiary mass that exceeds the secondary would consider- 
ably change the outcome of the evolution, because in that 
case it is the outer star that ascends the giant branch and 
possibly fills its Roche lobe before the inner binary has 
turned into a LMXB. The evolutio n of such triples has 
not yet been subject to any study (jEggleton fc Kiseleval . 
1199ft lEggletonll2006[ ). 

When the massive primary of the initial trip le ascends 
the red-giant branch a common-envelope ensues dWebbink , 
1984 iTaam fc Sandauisti l2000t iDelove fe Taaml . I201C ; 



Ge et a/J . l20ioir in which the inner secondary and the de 



generate core of the giant spiral in towards each other. 
A few Myr later, the stellar core explodes as a super- 
nova, forming a neutron star with a low-mass compan- 
ion in an eccentr ic orbit with a period of a few to sev- 
eral tens of days ([Delove fc Taaml [20101. The outer orbit 
is likely influenced both by the mass shell ejected in the 
common envelope, as well as by the supernova explosion. 
The mass shell containing 6 to 10 M©, may slow down 
the outer companion star, which may cause its orbit to 
shrink, but the removal of mass from the system may also 
widen it. The supernova explosion causes the inner binary 
to rec eive a velocity kick (jBlaauwl 119611 : iDewev fc Cordesl . 
Il987t ) which changes the orbits and may disrupt the triple. 
The kick results from two subsequent effects; the mass loss 
in the explosion and the intrinsic velocity kick imparted 
to the newly formed neutron star (see § 15.21 for details) . 
In the case of an electron-capture supernova the latter ef- 
fect is expe cted to resu l t in a relatively low-velocity kick 
< 50km/s (|Dewi et all 120051 ). which generally sufficices 
to keep the outer orbit bound by compensating the ef- 
fect of the mass loss in the supernova. In the surviving 
systems mass transfer will ensue in the inner binary af- 
ter several Gyr. This rapidly circularizes the inner orbit 
by tidal forces and turns the system in a Low Mass X- 
ray Bina ry (LMXB) resembling the bright Galactic bulge 
LMXBs ([Webbink et a/.IJl983t iTaaml . 119831 ). An example 
of such a triple is 4U 2129+47 (V1727 Cyg), in which a 5.24 
hour LMXB is accompanied by a spec tral type F-dwarf in 
an eccentric orbit of about 175 days ([Garcia et all Il989t 
iBothwell etall , 120081: iLin et all [20091) . For our model we 
require that, contrary to 4U 2129+47, at the moment of 
RLOF the LMXB has a period of at least ~ 1 day and 
will evolve to longer periods, rather than shrink due to 
the emissio n of gravitational waves ()Pvlvser fc Savoniid . 
119881 11989D . In the case of 4U 2129+47, the inner orbit 
will not expand due to the mass transfer and therefore will 
not perturb the outer orbit. 



Mass transfer in the inner binary causes the accreting 
neutron star to be spun up to a millisecond spin perio d 
(jAlpar et all 119821: iRadhakrishnan fc Srinivasanl . I1982D . 
After several tens of millions of years, the expansion of 
the inner orbit leads to a dynamically unstable situation 
with respect to the orbit of the outer star. Depending 
on the precise configuration at the moment when the sys- 
tem becomes dynamically unstable, either the donor of 
the LMXB, or the outer tertiary star or the neutron star 
can be ejected from the unstable triple. In the latter case 
a single MSP results, in the first-mentioned case a wide 
eccentric millisecond binary pulsar with a G-dwarf com- 
panion is produced, resembling J1903 + 0327. 

3. CONSTRAINING THE BIRTH CONDITIONS FOR 
J1903 + 0327 

We investigate the conditions under which the pro- 
posed triple scenario produces a system similar to 
J1903 + 0327. With the currently observed orbital pa- 
rameters of J1903 + 0327 (a ~ 1 23.2 R^, e ~ 44, m G ~ 
1.03 M© and M MSP ~ 1.67M© (jFreire et all I2010D ) the 
binding energy of the binary is about ~ 2.65 x 10 46 erg. 
This poses a minimum to the binding energy of the triple 
just before it became dynamically unstable, and allows us 
to calculate the orbital separation of both the inner and 
the outer orbit at the moment that the dynamical insta- 
bility sets in. In practice the binding energy of the stable 
triple will be higher than Ej, by the ratio of the masses 
of the e j ected star with respect to the triple, or ~ 20% 
([Heggid . 119751 : iFregeau et all I2004D . because the escaping 
star carries off some fraction of the binding energy to in- 
finity. 

The requirement that the binding energy of the triple 
must exceed Eb allows us to calculate the orbital separa- 
tion of both the inner and the outer orbit at the moment 
that the dynamical instability sets in. The criterion for 
dynamical stability sensitively depends on the separation 
and the eccentricity of the in ner and outer orbits, and o n 
the masses of the three stars (jMardling fc Aarsethl . |2001| ) . 
However, for the observed system J1903 + 0327 only the 
masses of two of the stars (the MSP and its current G- 
dwarf companion) arc known, as the initial secondary star 
was ejected from the binary. 

The initial donor in the LMXB must have been more 
massive than the current MSP companion in order to 
evolve first, but not so massive that mass transfer in the 
LMXB would be unstable; we therefore adopt a mass of the 
original close companion of the neutron star of 1.0 - 2.0 M© 
(In principle a secondary mass as low as 0.8 M© would suf- 
fice to warrant the inner binary to evolve into a LMXB, but 
in that case the tertiary star should be < 0.8 M©, which is 
smaller than than the observed 1.03 M©). For the neutron 
star we a dopt a mass o f 1.30 M^ before it starts accret- 
ing mass ([Schwab et all I2010D . and in the same range of 
masses, but less massive than the inner secondary star for 
the outer companion. In the specific case of J1903 + 0327 
the observed star is an 1.03 M© G-dwarf, which in our sce- 
nario would be identical to the initial tertiary star. The 
best match to the orbital characteristics are then obtained 
when we adopt the initial secondary mass to be ~ 1.1 Mq. 

During the LMXB phase mass is transferred from the 
inner companion (donor) to the neutron star. The time av- 
eraged mass-transfer rate can be estimated from the em- 
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pirical r elation based on t he initial orbital period of the 
LMXB dShore et oU 119941. In our model we assume that 
the neutron star accretes at most at the Eddington limit; 
for which we adopted m = 1.5 X 10~ 8 M©/yr. The left- 
over mass that is provided by the donor but that is not 
accreted by the neutron star is assumed to leave the in- 
ner binary with t he specific orbital angul ar momentum of 
the neutron star (Portcgics Zwart, 1995]). The currently 
observed mass of the neutron star is 1.67 M©, which in- 
dicates that in our model it must have accreted about 
0.37 M©, and that the donor star must have lost at least 
that same amount of mass. Mass transfer in such a binary 
system causes the orbit to expand quite dramatically, in 
particular if the mass transfer was not conservative. 

We can now determine the separation of the inner and 
the outer orbit for the initial triple (before the LMXB- 
phase) by an iterative procedure under the constraints 
that the total binding energy of the triple at the moment 
it became dynamically unstable is known, and by taking 
the effect of the mass-transfer process during the LMXB 
phase into account. Here we assumed that the mass lost 
from the LMXB leaves the triple in the form of a stel- 
lar wind. The resulting orbital constraints of this iter- 
ative procedure are presented in Fig. [I] as the solid and 
dashed curves; these give the most likely initial conditions 
for forming J1903 + 0327. 

With the semi-analytic procedure just described we still 
cannot predict the consequence of the dynamical instabil- 
ity, in particular because the final stage of the triple is 
highly dynamic and the orbital parameters do not allow 
us to predict the identity of the ejected star from first 
principles. However, one can imagine that it is quite likely 
that the less massive donor (< 0.73 M© when we adopt 
an initial secondary mass of 1.1 M©) is ejected, rather 
than the more massive MSP (~ 1.67 M©) or the outer 
star (~ 1.03 M©). At the moment when the dynamical 
instability sets in the donor still has part of its Hydrogen 
envelope, which explains why its mass exceeds that of the 
degenerate Helium core. 

4. SIMULATING MASS TRANSFER IN A TRIPLE SYSTEM 

We quantify our proposed formation channel for 
J1903 + 0327 by performing extensive computer simu- 
lations of triple systems, starting at the onset of the 
LMXB phase. The calculations ware performed us- 
ing the Ast rophysics Multipurpose Softw are Environment 
fAMUSEfl(P ortegies Zwart et aLl . l200l . 

In the AMUSE environment we resolve the dynami- 
cal evolution of the triple-star by a specialized numeri- 
cal orbit integration, which is written in FORTRAN77. 
The mass transfer in the inner binary is implemented in 
python. The outer star is treated as a point-mass and 
was not evolved during the simulation. The coupling be- 
tween the numerical orbit integration of the three stars and 
the stellar evolution ca l culati ons is realized using AMUSE 
(|Portegies Zwart et aZ.I . [2009h . The most important role of 
AMUSE is converting the units and to realize the commu- 
nication between the two codes. The former is done with 
a specialized unit conversion module and for the latter we 
spawned the different processes using the Message Passing 
Interface (MPI) to transfer the required data in the proper 

a see http : / /www . amusecode . org 



units and in discrete instances between the two codes. 

4.1. Numerical Method 

The orbits of the triple stars were integrated us- 
ing a regularized versio n of the Burlish-Stoer integrator 
(|Aarseth fc Zard . Il974al fbh , keeping the numerical error at 
machine precision, and allowing a maximum relative en- 
ergy error of O(10~ 14 ) per integration of the outer orbit. 

During the dynamical evolution we resolved the mass 
transfer and consequential change in orbital paramters of 
the LMXB that is orbited by the outer star. 

The rate of mass-transfer in the inner binary is cal- 
culated from the empirical relation based on the orbital 
period at the onset of RLOF in the LMXB: (rh) = 6 x 
10 - 10 (P OT b (initial) / lday) M© /yr dShore et all fl99l . 

The neutron star was allowed to accrete at most at the 
Eddington rate and any surplus mass is assumed to leave 
the binary with the angular momentum of the accreting 
neutron star, but lost adiabatically from the triple. 

We perform the mass transfer in the inner binary every 
time the outer orbit has had 10 revolutions, after which the 
numerical orbit integration was updated using the newly 
calculated orbital parameters which resulted from carry- 
ing out the mass transfer. We varied the interval between 
which mass transfer in the inner binary was conducted 
between every 1 to 1000 outer orbits, but this choice did 
not significantly affect the results. (In S I5. 31 we present 
the results of a series of simulations where we decoupled 
the mass-transfer process from the gravitational evolution 
after the triple has become dynamically unstable.) 

4.2. Results of the simulations 

We initialized 10 3 binaries and calculated the evolu- 
tion for each up to an age of at most 10 Gyr or until 
the mass of the donor star drops to the mass of the de- 
generate helium core for a population-II star ( £ 0.4 M©, 
iTauris fc Savoniiel I1999D . 

Each triple was initialized by randomly selecting the ec- 
centricity of the outer orbit from the thermal distribution 
and the initial separations of the inner and outer orbits for 
each triple-evolution calculation are selected using the it- 
erative procedure described in §[3] (see the thick solid and 
dashed curves in Fig.[T]). For each simulated binary we 
randomly selected the inclination of the inner orbit with 
respect to the outer orbit, the longitude of the ascend- 
ing node, the argument of periastron and the phases of 
the two orbits. The preference in the orbital elements in- 
troduced by the supernova explosion may affect the rate 
of ejected MSPs relative to those that stay in a binary, 
but we ignore that complication. We found no significant 
correlations between the final outcome of the simulations 
and longitude of the ascending node, the argument of pe- 
riastron or the phases of the two orbits. The anisotropic 
velocity caused by the mass loss in the supernova of the 
inner binary therefore is not expected to have a significant 
effect on the surviveability of the triple. 

During the evolution of the triples, small eccentricities 
( ^ 0.1) are commonly induced in the inner orbit, in par- 
ticular when the triple approaches the regime where it be- 
comes dynamically unstable. We stop a simulation when 
the eccentricity of the inner (LMXB) orbit exceeds 0.3 for 



4 



more than 10 5 years. Such high eccentricities can be in- 
duced shortly by a strong interaction with the outer star, 
which typically results in the break-up of the triple, or 
by a long term se cular p e rturb ations of the inner orbit 
by the outer star (|Kozail . I1979D . In the latter case or- 
bital variations result naturally from the secular evolution 
of the triple and do not directly lead to a dynamical in- 
stability We still decided to terminate such simulations 
because it becomes hard to follow the mass-transfer pro- 
cess within the inner binary. It would require extensive 
hydro-dynamical simulations to study the consequences 
of a Kozai resonances in a triple with a Roche-lobe fill- 
ing inner binary. About 10% of our simulations ware 
stopped as a consequence of this effect. Mass transfer 
in eccentric orbits is generally ill understood, although 
courageous attempts a re un derway to tack l e this p roblem 
(ISepinskv et all l2009t I20I0L lLaioie fc Siilsl . l2010albh . 

The majority (917) of the binaries become dynamically 
unstable long before the other stopping criteria apply, in 
which cases we continue to resolve the dynamics by inte- 
grating the equations of motion and resolving the internal 
mass transfer until one of the stars escapes. The orbital 
separation of the LMXB at which the triple is expected to 
become unstable is indicated by the dotted curve in Fig. [I] 

An illustrative example of the evolution of the period 
of the inner and outer orbits is presented in Fig.O The 
evolution of the inner LMXB also drives the expansion of 
the outer orbit, in particular by the dynamical coupling 
between both orbits and in a lesser extend by the mass 
lost from the inner LMXB in those cases that mass trans- 
fer proceeds non-conservatively. This relative softening 
causes the final MSP binaries to be somewhat wider than 
observed in J1903 + 0327, and as a consequence the triple 
remains stable for somewhat longer than expected based 
on our analytic energy balance (see §[3j) , and eventually re- 
sults in the MSP to be somewhat more massive (by about 
0.2 M Q ) than observed in J1903 + 0327. We can compen- 
sate for this by adopting a slightly (~ 20 %) smaller outer 
orbital separation at the onset of the LMXB phase. 

The resulting parameters of the MSP with G-star bi- 
nary are presented in Fig.[3j and straddle the observed 
orbital separation and eccentricity of J1903 + 0327. The 
phase of mass transfer for these binaries lasted for about 
23.4 ± 9.3 Myr, after which the donor was ejected in about 
one-fifth ~ 20% of the cases. The resulting binaries, for 
those with a < 10 4 R©, had an average orbital separation 
of 175 ± 145 Rq, and an eccentricity of 0.65 ± 0.19. Dur- 
ing the LMXB phase, the neutron stars were able to grow 
from 1.30 M© to 1.64 ± 0.12 M Q . The mean mass of the 
ejected donor was 0.72 ± 0.13 M©. 

Our success in reproducing the observed parameters of 
the binary MSP J1903 + 0327 demonstrates that its pro- 
genitor may well have been born as a triple star. How- 
ever, the here described pin-pointed search of parame- 
ter space makes it impossible to derive the birthrate for 
J1903 + 0327, which we will calculate in the next §. 

5. HOW MANY J1903 + 0327-LIKE SYSTEMS ARE THERE IN 
THE GALAXY? 

We calculate the formation rate of J1903 + 0327-like 
systems by determining their birth rate with respect to 
that of ordinary LMXBs. The reason for this approach 
is the great uncertainty in the number of LMXB progen- 



itors because binaries with such extreme mass ratios are 
observationally unknown. 

This ratio depends on several factors 

1 . The number of triple systems with suitable parame- 
ters. 

2. The fraction of triples that ensues and survive the 
common-envelope phase and the supernova explo- 
sion in which the neutron stars is formed 

3. The fraction of those triples that lead to 
J1903 + 0327-like systems (instead of single MSPs, 
classic MSP binaries or MSP binaries with a triple 
companion) 

We will discuss each these factors below, in § 15.11 § 15.21 
and § 15.31 respectively. 

5.1. The number of suitable triple stars 

According to the Hipparcos database the ratio of hier- 
archical higher order multiple stellar systems (404) to bi- 
naries (1438) is 404/1842 ~ 22 and the outer star have a 
rathe r flat mass distribution (Eggl eton fc Tokovininl [2008, 
We assume that for (inner) binaries that are pro- 
genitors to LMXBs the same ratio holds, even though there 
are none in this catalogue. We further require that the 
mass of the outer stars is about ten times lower than that 
of the inner binary. We then find a ratio of suitable triples 
to LMXB progenitor binaries of a few per cent. We present 
in § 16.21 observational evidence that indeed such triples ex- 
ist. 

5.2. The fraction of triples that survive the supernova 

explosion 

We will not dwell on the details of the common enve- 
lope evolution, but assume that it leads to the spiral-in 
of the inner two components without much affecting the 
outer star. The consequences of the common envelope do 
not qualitatively affect our result, but have a strong effect 
on the derived birth rate (see § 15. 4p . In § 16.21 we discuss 
that the orbit of the tertiary star that orbits the LMXB 
4U 2129+47 may have experienced a considerable reduc- 
tion during the common-envelope, but we consider this in- 
sufficient evidence to draw general conclusions regarding 
the effect of the common envelope on the outer orbit. 

The effect of the supernova explosion in the inner binary 
on the orbital parameters of the triple can be calculated 
relatively straight-forward, and has a profound effect on 
the survivability of the triple because the weakly bound 
outer orbit is easily disrupted. 

The mass loss in the su pernova explos i on causes the in - 
ner binary to be ejected (jBlaauwl 119611 : iBoersmal Il96ll ). 
In order to keep the triple bound a small asymmetric ve- 
locity kick i mpar t ed to the newly formed neutron star 
is required (iHillsl 119831 : iBrandt fc Podsiadlowskl 119951: 
iTauris fc Takensll 1998ft . In Fig.g] we show the survival 
probability of several combinations of inner and outer or- 
bits as function of the asymmetric kick magnitude. In 
these calculations the effect of the Blaauw-Boersma kick 
was self-consistently taken into account to calculate the 
survival probability. This fraction ranges from ~ 40% for 
small kicks to zero for kicks significantly above 100 km s _1 . 
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To study the future evolution of the triple through the 
LMXB phase, we need to derive the orbital parameters 
as a function of the most probable kick velocity and the 
amount of mass lost in the supernova. We consider several 
combinations of both, and the distinction between them 
can be made by the mass of the initial primary star and 
the moment it filled its Roche-lobe. 

In §[3] we discussed the most likely range of parameters 
at birth, and we adopt these same numbers here to study 
the surviveability of the triple. Here we make a distinction 
between two types of supernovae; in one case we adopt 
that the inner primary at the moment of the supernova 
is a 1.9 Mq Helium star with a degenerate ONeMg core 
that is the left-over of the initial 9-12 M Q primary star. If 
the initial primary star was born somewhat more massive, 
10-13 M Q and was stripped from its Hydrogen envelope 
at a later stage of its evolution it's Helium core may have 
grown to ~ 2.7 Mq. The ONeMg star that experiences 
an electron-capture supernova loses ~ 0.6 Mq and if the 
exploding star has a more massive Helium core the mass 
loss is ~ 1.4M Q . In both cases we adopt a neutron star 
mass of 1.30 Mq. 

Apart from the smaller mass lost in the explosion 
of the ONeMg core, the neutron star is also expected 
to receive a smaller velocity kick upon birth. For 
these electron capture supernovae we adopted a Gaus- 
sian distributio n for the velocity kick with a disper- 
sion of 20 km/s (|Podsiadlowski et all 12004 IScheck et all 
12004 in a random direction. From studies of single ra- 
dio pulsars in the solar neighborhood s everal kick veloc- 
ity d i stributions have been con s tructed dLvne fc Lorimerj . 
1994lHansen fc Phinnevl 119971: ICordes fc Chernofil . 11998: 



Arzoumanian et cdl I2002T ) , most of them with considerbly 



higher velocity than for the ONeMg supernovae. These 
latter kick velocity distributions are expected to be more 
suitable for neutron stars formed from an isolated star or 
a non-interacting binary, in which case t he neutron star is 
formen by the collapse of a n iron core (|van den Heuvell . 
12004 iPodsiadlowski et all 12004 . These kicks range 
from single (zero-centered) Gaus sian distribut i ons with 
a velocity dispersion of 265 km/s iHobbs et al\ (|2005l ) to 
more complicated distributions like the proposed two 
Gaussians with dispersions of 175 km/s and 700 km/s 
and relative probabilty o f 0.86 and 0.14, respectively 
(|Cordes fe Chernoffl . 119981 ). 

We study the probability that a triple survives the su- 
pernova by means of Monte-Carlo simulations. The sim- 
ulated triples had the following characteristics: The inner 
binary consists of a 1.9 Mq Helium star with a degener- 
ate ONeMg core or a ~ 2.7 Mq Helium star, and a 0.8 - 
2.0 Mq secondary. The latter was selected randomly with 
equal probability within the interval. The orbital sepa- 
ration of the circular inner binary was taken flat in log 
between IRq and 100 Rq. We adopted these parame- 
ters from the populatio n synthesis calculation of LMXBs 
l|Willems fc K olb. 2003j), in parti cular using the results o f 
their models KM25 to KM100 of (jWillems fc Kolbl . I2003D . 

The mass of the tertiary star was selected to be less mas- 
sive than the secondary, but not smaller than 0.8 M Q (in 
theory there is no lower limit for the mass of the tertiary 
star). The outer orbital separation was chosen with a prob- 
ability distribution flat in log with a maximum of 10 4 Rq. 
The minimum separation was chosen to be consistent with 



a dynamically stable initial triple (adopting a 10 Mq pri- 
mary and an inner orbital separation of 200 Rq) and the 
earlier selected secondary and tertiary masses. The lat- 
ter two stars ware assumed not to accrete any material 
throughout the common envelope and supernova explo- 
sion. The eccentricity of the outer orbit before the super- 
nova was selected at random from the thermal distribution 
between a circular orbit and a maximum which was cho- 
sen such that the triple is dynamical stable. The selection 
of the minimum orbital separation of the outer star and 
it's eccentricity therewith becomes an iterative procedure. 
The other orbital elements are selected randomly, as we 
described in § 14.21 

For each system we calculate the effect of the com- 
bined Blaauw-Boersma and intrinsic velocity kick on 
the inner and the outer orbit, the latter kick was as- 
sumed to be isotropic. In our simulations we varied 
the mass loss in the supernova and the velocity dis- 
tribution of the asymmetric kick . For electron-capture 
supernovae (|Podsiadlowski etall 12004 the fraction of 
surviving triples is 1/3. For higher kick velocities we 
adopted that the exploding star was 2.7 M and as 
a consequence the fraction of s urvivors drop s to 1/25 
(I Arzoumanian et aZ.Ll2002tj. 1/2 8 (|Hobbs et a/.l . l2005l) and 
1/50 (|Cordes fc Chernofil . 119981) where the quoted litera- 



ture refers to the adopted kick velocity distribution. Note 
here that the smaller amount of mass lost in the electron- 
capture supernova explosions helps considerably in pre- 
serving more triples, as opposed to the more violent kicks. 

We conclude that for every 3 inner binaries that sur- 
vive the electron-capture supernova the outer tertiary star 
remains in orbit around the inner binary, but that this 
fraction may drop considerably (to 1/50) when more mass 
is ejected in the supernova shell and the kick velocity is 
higher. Varying the amount of mass lost in the super- 
nova explosion has a profound effect on the survivability 
of the triple, in particular since the Blaauw-Boersma kick 
imparted on the inner binary is proportional to this mass 
loss, and to the relative orbital velocity of the inner binary. 
It is interesting to note that the vast majority of the triples 
that survive the supernova explosion are dynamically sta- 
ble, but their orbital eccentricity tend to be considerably 
higher than according to the thermal distribution. 

5.3. The fraction of surviving triples that lead to systems 
like J1903 + 0327 

We synthesize the Galactic population of binaries like 
J1903 + 0327 by randomly selecting 10 4 triples that sur- 
vived the supernova explosion of the previous §, and con- 
tinue their evolution in AMUSE (see § 14. 1[) . Instead of 
performing a self-consistent evolution as adopted in § 14.21 
to validate the proposed scenario, we tentatively decou- 
pled the mass transfer process from the orbit integration. 
The population synthesis simulations start by resolving 
the mass transfer in the inne r binary until the triple be - 
comes dynamically unstable ([Mardling fc Aarsethl . I200l1 ). 
after which we continue the simulation by resolving the 
dynamics of the 3-body system until one star is ejected. 
During this latter part we ignore the mass transfer. 

Note that this decoupled approach, though computa- 
tionally cheaper by about a factor of 10 , is not a-priori 
less reliable than the self-consistent simulations in § 14.21 
because our numerical methods provide no self consistent 
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way to resolve the mass transfer process in eccentric orbits. 

The majority ( £ 90%) of our simulations lead to the 
disruption of the triple, which results in the ejection of the 
initial primary (MSP), the secondary (partially stripped 
sub-giant or white dwarf) or the tertiary (main-sequence 
F-, G- or K-dwarf) star, leaving the other two stars in a bi- 
nary. The respective ratio at which these occur are 0.013, 
0.683 and 0.304 for circular outer orbits (which is unlikely 
after the supernova) to 0.03, 0.489, and 0.481 for highly 
eccentric (e ~ 0.9) orbits. For extremely high eccentrici- 
ties (e ~ 0.99) these fractions become 0.05, 0.75 and 0.20, 
respectively. The outer orbits tend to have high eccen- 
tricities because they often barely survive the supernova 
kick, and as a consequence the total fraction of systems 
in which the MSP remains in a binary with the original 
outer component is £ 0.5; the fraction of single MSPs is 
- 0.05. 

The two methods through which we resolved the triples 
(self-consistent, see § 14.11 as opposed the here adopted de- 
coupled approach) give slightly different results, but these 
can be related to the variation of the implementations. 
The biggest effect is illustrated in Fig.[5J where we demon- 
strated how the slow changes in the inner orbit drive the 
secular evolution of the outer orbit. The main conse- 
quences of this coupling are the higher mass of the neutron 
star by the time the triple becomes dynamically unstable 
and the small fraction of triples that survive the entire evo- 
lution because the adiabatic expansion of the outer orbit 
prevents the triple from becoming dynamically unstable. 
The consequential loss of systems however, is well compen- 
sated by the increased probability that the donor in the 
LMXB is ejected when the triple becomes dynamically un- 
stable. The differences between the two numerical imple- 
mentations affect our estimates for the birthrate (see § 15. 4[) 
on a £ 20% level. For clarity, the statistical uncertainties 
between the two different numerical approaches depend on 
Poissonian arguments rather than on the lack of our un- 
derstanding of parts of the physical process, such as the 
common envelop evolution, the process of mass transfer 
in non-circular orbits and the non-linear effects in the dy- 
namically unstable configuration just before the triple is 
resolved. We control the latter by adopting a high order 
and numerically extremely precise algorithm (see § 14. ip . 
However, the tipping-point physics of ejecting the MSP, its 
close white-dwarf companion or the main-sequence outer 
star remains elusive, even when integrating near machine 
precision. 

5.4. The number of J1903 + 0327-like systems in the 

Galaxy 

The ratio of the birth rates of ordinary LMXBs to sys- 
tems like J1903 + 0327 is a combination of the factors 
derived above: a few per cent for the number of suitable 
triples, a reduction of a factor of ~ 3 owing to the effect of 
an electron-capture supernova, and then a fraction ;$ 0.5 
of systems in which the inner secondary is ejected. For 
larger average kick veloc ities ([Cordes fc Chernofl [l998; 
lArzoumanian etHHl , l2002f ) the fraction of triples that sur- 
vive the supernova drops from 1/3 to 1/50 (see § 15. 2|) . The 
higher kicks cause an even stronger reduction in the num- 
ber of single MSPs, because outer orbits tend to be highly 
eccentric, which leads to a reduction of the probability that 
the MSP is ejected once the triple becomes dynamically 



unstable. The birth rate for systems like J1903 + 0327 
is then in the range of 2 x 10~ 4 to 3 x 10~ 3 times the 
birth rate of ordinary LMXBs. Estimates of the total num- 
ber of LMXBs in the Galaxy are in the range 10 4 - 10 5 
(jCote k Pvlvseii 119891) . The typical life time of a LMXB 
is of the order of a Gyr, whereas the life time of systems 
like J1903 + 0327 is limited by the life time of the MSP 
and the difference between that of the inner secondary star 
and the outer tertiary star. Both stars have a mass of the 
order of 1 M Q and their lifetimes exceed several Gyr, which 
is a factor of a few longer than the lifetime of the LMXB 
phase. The number of J1903 + 0327-like systems in the 
Galaxy then is at least 30-300 for electron-capture super- 
novae but drops to 3-30 when we adopt considerably more 
mass to be lost in the supernova and higher ve l ocity kicks 
(jCordes fc Chernofl 119981: lArzoumanian et "oD. 12002ft . 

6. DISCUSSION 

We went through considerable effort to explain the ex- 
istence of J1903 + 0327 as the result of the complex evo- 
lution in a hierarchical triple star system. In this section 
we will argue that though, unlikely as it sounds, our pro- 
posed model is currently the only viable model available, 
and that all other existing models fail to explain all the 
characteristics of J1903 + 0327. After that we will indi- 
cate a number of other sources that have evolved in quite 
a similar way, and that therefore support the triple sce- 
nario. 

6.1. Shortcomings of earlier models for J1903 + 0327 

In the tripl e scenario for the f ormat ion of J1903 + 0327 
proposed by IChampion et all (|2008ft . the observed ~ 
95 day orbital period is that of an unobserved massive 0.9- 
1.1M white dwarf, while the observed G-dwarf has a con- 
siderably wider orbit. The observed eccentricity of the in- 
ner orbit would in this case be driven by th e secular evolu- 
tion via the Kozai mechanism (jKozail Il962ft . This scenario 
has a number of serious shortcomings, one of which is the 
unusually high white-dwarf mass, which leaves little room 
to produce a massive 1.67 M Q neutron star. Since it's dis- 
covery, the observed change in eccentricity e ~ 10~ 16 s -1 is 
three orders of magnitud e smaller than prediction s based 
on the Kozai mechanism (jGopakumar et all , 120091 ). which 
excludes the presence of a highly inclined tertiary star. We 
exclude this scenario therefore from further consideration. 

An alternative to the previous model would be the di- 
rect formation of a MSP in a supernova explosion, or 
via t he fall-ba c k of material in a circum neutron-star 
disk (|Liu fc Lil . I2009D . A variety of arguments against 
the direct formation of rapidly spinning pulsar with a 
sm all surface magnetic fi eld (2 x 10 8 G) was provided 
by IChampion et all (|2008ft . These include: first, of the 
50 neutron stars in young supe rnova remnants, none ar e 
fast-spinning low- field pulsars (jKaspi fc HelfancH . 120021 ). 
Second, a "born- fast" scenario for J1903 + 0327 would 
likely account for the 15 isolated MSPs detected in the 
galactic disk, but the spin distribution, space velocities 
and energetics of these single MSPs are indistinguish- 
able from those of recy cled, not "born fast" binary MSPs 
([Archibald et all , 12009ft : while their space velocities and 
scale he ights do not match those of non-recycled single 
pulsars ([Lorimer et all , 12007ft . Third, magnetic fields in 
young pulsars likely originate either from dynamo action 
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in the proto-NS (jThompson fc Duncanl . I1993T ) or through 
compression of "frozen-in" fields of the progenitor star dur- 
ing collapse ( Ferrario fc Wickramasinghel 1200 6) . and no 
young pulsars with B < 10 10 Gauss, like J1903 + 0327, 
are known. We add a fourth argument: before the direct 
collapse, the NS progenitor ascends the giant branch star 
and common-envelope evolution with the 95 day-orbit G- 
dwarf would have dramatically reduced the orbital period, 
in the same way as normal LMXBs form. For these four 
reasons, a core-collapse "born-fast" MSP can be ruled out. 

The above arguments hold equally well against the 
accretion-induced collapse of a massive and rapidly ro- 
tating white dwarf to a neutron star (NS). While this 
may lead to millisecond rotation periods, the NS formed 
here will be just as hot and differentially rotating during 
the early liquid phases as a NS formed by core collapse. 
In both cases 10 53 erg in gravitational energy is released, 
which will erase all memory of the violent mechanism in 
which the NS was formed. Therefore, dynamo action in 
both cases will not depend on the formation mechanism, 
and there is no fundamental reason for expecting a weak 
magnetic field in NSs formed by accretion-induced col- 
lapse. Furthermore, the direct accretion-induced collapse 
of a white dwarf to a MSP in population studies produces 
binaries with a orbital period 20 days, whic h is consider- 
ably shorter than observed in J1903 + 0327 (|Chen et all 
2010). In addition this scenario requires the white dwarf to 
accrete from an evolved companion star, which is inconsis- 
tent with the observed main-sequence companion G-star. 
This model thus fails to explain J1903 + 0327. 

Finally, a formation scenario where J1903 + 0327 is 
formed when the donor star in the inner binary is ab- 
lated and destroyed, like is the case for the ''black- 
widow" system PSR B1957+20 (jFruchter et all Il988f) . 
has thee problems. First, the observed timescales for 
straightforward evapora tion of the donor star are too long 
(jChampion et aLl , l2008t h Second, formation o f such a sys- 
tem likely involves an exchange interaction (|King et all 
l2003h . which would be greatly impeded by the outer G- 
dwarf companion. Third, even if no exchange took place, 
the slow evaporation of the donor in the inner binary re- 
quires the triple to be dynamically stable with respect to 
the G-star in its current orbit. Since the supernova ex- 
plosion can at most account for a reduction of a factor of 
2 in the orbital separation, the common-envelope should 
in that case be responsible for a further reduction from 
the initial orbital separation of XL 560 Rq to the currently 
observed ~ 123 Rq, which requires considerable fine tun- 
ing. We therefore conclude that J1903 + 0327 is unlikely 
to have originated through a black-widow-like scenario. 

Each of the above models has serious shortcomings, 
and we conclude that none of the scenarios discussed 
above give a satisfactory explanation for the formation of 
J1903 + 0327. 

6.2. The missing link 

There are currently no triples known with parameters 
suitable for evolving into systems like J1903 + 0327. How- 
ever, there is observational evidence that such triples ex- 
ist, as it is proposed that 4U 2129+47 (V1727 Cyg), a 
5.24 hour LMXB, is accompanied by a sp ectral type F- 
dwar f in an eccentric orbit o f about 175 days dGarcia et all , 
Il989t iBothwell et all l2008t iLin et all |2009| ). This triple 



can have formed in the same way as J1903 + 0327 with 
the exception that after the common-envelope and the 
subsequent supernova explosion the inner binary period 
was smaller than the b i furcat i on pe riod, of about one 
day dPvlvser fc Savoniid . 119881 11989ft . The consequence 
of such a short orbital period is that the LMXB evolves to 
an even shorter orbital period. Interestingly the 175 day 
orbit of the outer star is too small to have been dynami- 
cally stable at the birth of the triple ([Mardling fc Aarsethl 
I2001h . We argue that the common-envelope and/or the 
supernova may have reduced the orbital separation of the 
outer star. 

We validated the probability that a binary orbit shrinks 
as a result of the supernova by means of population syn- 
thesis and conclude that in a fraction of 0.511 of the 
triples that survive the supernova the separation of the 
circularized outer orbit is smaller than the initial orbit. 
However, if the progenitor of 4U 2129+47 had parame- 
ters comparable to what we derived for J1903 + 0327 in 
§[2J the common-envelope phase of the inner binary must 
have resulted in a reduction of the outer orbit as well, by 
^ 210 Rq. Such a reduction in the separation of the outer 
orbit by the common envelope enormously boosts the sur- 
viveability of the triple in the supernova explosion (see 
Fig. 2]) , and therefore dramatically increases the birthrate 
of binaries like J1903 + 0327. 

Following the same scenario as for forming 4U 2129+47 
but with a less massive ( ^ 8Mq) initial inner primary 
star the inner binary could evolve into a cataclysmic vari- 
able such as EC 19314-5915, with an orbital period of 
4.75 hours. The observed radial velocity of ~ 9km/s has 
been attribute d to a G8-dwarf in orbit around the CV 
(Bu ckley et all 17992 i. which is consistent with a semi- 
major axis of ~ 2400 R© 

The real missing link would be the discovery of a rela- 
tively wide ( ^ 30 R©) LMXB that is orbited by a tertiary 
low-mass main-sequence star. The mass-transfer phase in 
our simulations averaged about 23Myr (see § I4.2[) , while 
the lifetime of the MSP in J1903 + 0327 is at least 1 Gyr. 
We therefore expect that the Galaxy contains at most 7 
such wide triple LMXBs. 

7. CONCLUSIONS 

We discussed the evolution of triple star systems 
through a range of dramatic events, including several tidal 
circularizations, a common-envelope phase, a supernova 
and a stable phase of mass transfer that eventually leads 
to a dynamical instability in which one star is ejected. 

In particular for producing a system like J1903 + 0327 
we require a triple to be born as a rather ordinary dy- 
namically stable hierarchical system of which the inner 
binary consists of a 9-13 M Q and a 0.8-2.0 M star in a 

200 Rq separation. This binary is orbited by a main se- 
quence star with a mass smaller than the initial secondary 
(< 2.0M Q ) with a semi-major axis £ 560 Rq. The chance 
that the triple survives the inevitable chain of events is not 
large but the result is profound and provides a satisfactory 
explanation for a number of known systems in the Galaxy, 
including J1903 + 0327, 4U 2129+47 (see §EU and EC 
19314-5915. 

The range in possible observable stages in the evolution- 
ary sequence for forming a system like J1903 + 0327 sen- 
sitively depends on the orbital separation of the inner bi- 
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nary after the supernova, which makes the distinction be- 
tween evolving into a binary MSP like J1903 + 0327 or a 
LMXB with an outer tertiary companion like 4U 2129+47. 
The evolution of the triple naturally leads to a cataclysmic 
variable like EC 19314-4915, if in addition to a short post- 
common envelope period the initial primary star evolves 
into a white-dwarf rather than a neutron star. We expect 
that such triple CV's are rather common in the Galaxy. 

We demonstrated that our model can indeed reproduce 
J1903 + 0327 qualitatively, and we estimate the number 
in the Galaxy by performing extensive population synthe- 
sis of post-common envelope triple systems (see §[5j - Our 
starting conditions are the progenitor of a LMXB which 
is accompanied by a third low-mass companion in a rel- 
atively wide orbit. Based on the observed statistics for 
such systems, their survival in the electron-capture super- 
nova explosion of an ONeMg star and the consequences 
of the dynamical instability which results from the mass 
transfer in the inner binary, we conclude that the for- 
mation rate of J1903 + 0327-like systems is ~ 3 x 10 -3 
times that of LMXBs. With a life time at least as long 
as that of LMXBs, and an est imated total number o f 10 4 
- 10 5 LMXBs in the Galaxy (|Cote fc Pvlvseii Il989ft . we 
expect at least 30-300 systems like J1903 + 0327 in the 
Galaxy and an order of magnitude smaller number of single 
MSPs. The longer lifetime of the MSP binary compared 
to LMXB's results in an increase of this number of a factor 
of a few. In the most pessimistic scenario, when we adopt 
a higher velocity kick, this number drops to about 3-30 
MSP binaries like J1903 -I- 0327 in the Galaxy, and a few 
single MSPs. 



With a birthrate for Galactic LMXBs of 3.2 x 
!Q- 6 /vr (iKalogera fe Webbinkl . 119981) to 7 x 10- 6 /yr 
(|Cote fe Pvlvserl 1989D we conclude that systems like 
J1903 + 0327 form at a rate of < 2.1 x 10" 8 , and a ten 
times smaller rate for single MSPs. These low rates makes 
the proposed scenario unlikely, as we already expected, 
but sufficiently probable that the Galaxy should contain a 
few tens to hundreds of objects with characteristics similar 
to J1903 + 0327 and consequently provides a satisfactory 
explanation for J1903 + 0327. Since the birthrate of single 
MSPs is expect ed to be quite similar to that of LMXBs 
(|Dai fc Li I2010D our proposed triple scenario does not sig- 
nificantly contribute to the formation of single MSPs. 
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Fig. 1. — A triple with a 1.3 Mq neutron star and 1.1 Mq G-dwarf in a circular inner orbit (thick solid curve [red]) and a 1.0 Mq G-dwarf 
in an eccentric outer orbit (dashes) is stable and has the same binding energy as J1903 + 0327. These curves are calculated using the iterative 
procedure described in §[3] The orbit of the LMXB expands as a consequence of the mass that is transf erred from the initial secondary 
donor star to the neutron star. According to the analytic expression provided by (Mardling & Aarscth, 2001, their Eq. 90), the triple becomes 
dynamically unstable as soon as the orbital separation of the LMXB reaches the dotted curve. At this moment the triple dissolves and one of 
the three stars is ejected, leaving behind a binary with an orbital separation roughly somewhere between the dotted and the dashed (black) 
curve (see Fig.(3j- 
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Fig. 2. — Evolution of the orbital period of a rather typical case which evolves to a system very similar to J1903 + 0327. The square (red) 
represents the final orbit of this simulation, which in this particular case resulted in a binary consisting of a 1.73 Mq MSP and a 1.0 Mq 
companion in a ~ 100 day orbit with e ~ 0.43. The bottom line (red) indicates, as a function of time the evolution of the inner LMXB, 
whereas the top line (blue) represents the evolution of the outer orbit. 
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Fig. 3. — The eccentricity (e) and semi-major axis (a) of the MSP binaries with G-dwarf companions that resulted our from simulations 
(dots [red], see §[3j) . These binaries are initialized to mimic J1903 + 0327, and which are described in §[3] and represented by the thick solid 
and dashed curves in Fig.[T] and which are reproduced here to guide the eye regarding the most likely range of final binary parameters. The 
square (blue) indicates the current orbital separation and eccentricity of J1903 + 0327. Note that the final eccentricity of the surviving binary 
is unrelated to the eccentricity of the initial outer orbit (black dashed curve) . 
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Fig. 4. — Survival probability for a triple in which a 2.7 Mq star explodes in a supernova explosion to a 1.3 Mq neutron star. The inner 
companion star is 1.1 Mq and the outer star is 1.03 Mq. The inner orbit is circular and with a semi-major axis of 30 Rq (panel a) and a 
circular outer orbit is 1000 Rq, 3000, 9000 Rq (top to bottom). Panel b gives the supernova survival fraction for an outer separation of 
3000 Rq and an inner orbit of 10, 20, 90 Rq (bottom to top curve). Note that we adopted here high mass-loss in the supernova contrary 
to an electron-capture supernova, which suppresses the survival rate. 



